Dipole-Dipole coupled double Rydberg molecules 
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We show that the dipole-dipole interaction between two Rydberg atoms can give rise to long range 
molecules. The binding potential arises from two states that converge to different separated atom 
asymptotes. These states interact weakly at large distances, but start to repel each other strongly 
as the van der Waals interaction turns into a resonant dipole-dipole interaction with decreasing 
separation between the atoms. This mechanism leads to the formation of an attractive well for one 
of the potentials. If the two separated atom asymptotes come from the small Stark splitting of an 
atomic Rydberg level, which lifts the Zeeman degeneracy, the depth of the well and the location of 
its minimum are controlled by the external electric field. We discuss two different geometries that 
result in a localized and a donut shaped potential, respectively. 

PACS numbers: 78.60.-b,32.80.Rm,32.80.Rm 
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A gas of cold Rydberg atoms, atoms in states of high 
principal quantum number n, is more like a solid than a 
room temperature gas. On the 1 fis time scale of a typical 
experiment cold, 300/iK, Rydberg atoms at a density of 
10 9 cm -3 only move 2% of their average interatomic spac- 
ing. The atoms are effectively frozen in place. Equally 
important, the atoms interact, due to their large electric 
dipole transition moments fi n 2 , which lead to dipole- 
dipole interactions of magnitude Vdd = n 4 /i? 3 where R 
is the displacement from atom 1 to atom 2 We use 
atomic units unless specified otherwise. For n = 40 atoms 
at 10 9 cm~ 3 , Vdd = 3.8 MHz, and dipole-dipole energy 
transfer is easily observed in 1 fj,s [2j-|5|. It is closely re- 
lated to the motion of an exciton in a solid and to the 
Forster energy transfer observed in light harvesting sys- 
tems @, 0] ■ The dipole-dipole interactions of Rydberg 
atoms are also the basis of proposed approaches to con- 
structing quantum gates @, ■ Their second attraction 
is that it is possible to form novel long range molecules. 
In what are termed trilobite molecules, a Rydberg atom 
is bound to a ground state atom by the short range in- 
teraction between the Rydberg electron and the ground 
state atom [l^, 11 1. The trilobite molecules are huge, 
with an internuclear separation of n 2 . In addition, ev- 
idence for longer range molecular resonances has been 
presented, and stable macrodimers, formed by van der 
Waals interactions in weak electric fields, have been ob- 
served 12, 13j]. The macrodimers are huge, with internu- 
clear spacings exceeding 1 /im. However, they are fragile, 
existing only in very specific electric fields [13J . 

Here we demonstrate that it is possible to form robust 
long range molecules in which two Rydberg atoms are 
bound to each other by the dipole-dipole interaction. The 
essence of the idea is the following. If there are two sets 
of dipole-dipole potentials converging to two energetically 
close R = oo asymptotes, there are potentials converging 
to the upper asymptote which have wells due to repulsion 



from potentials converging to the lower asymptote. If the 
two R = oo asymptotes are produced b y li fting the Zee- 
man degeneracy using the Stark shift [14j | . the minima 
of the wells occur at extremely long range, R > 10n 2 . 
In other contexts the Zeeman degeneracy is an unwanted 
complication [lj| [l6|, but here it is an asset. A partic- 
ularly attractive feature of one of these wells is that, for 
any n, the location of its minimum depends on the Stark 
shift. It can, for example, be adjusted to match the sepa- 
ration of ground state atoms in a lattice [l?! , offering the 
prospect of a different way of making an ordered sample 
of Rydberg atoms [l8|,[l9|. 

After introducing the geometry of the problem we con- 
sider the simple case in which the two atoms are aligned 
in the field direction, R || E, for in this case the origin of 
the well is particularly clear. We then consider the more 
general case in which E and R are not parallel to show 
that the well has a finite extent in the angle 6 between 
E and R and that a well also exists for E _L R when 
the Stark shift is reversed in sign. Finally, we comment 
briefly on the properties, production, and detection of 
these molecules. 

The geometry of the two atom system under consider- 
ation is shown in Fig. [TJi. The electric field E is in the z 
direction, which we choose to be the axis of quantization, 
and R is at the angle 9 from the z axis. We here consider 
nsi/2 and np 3 / 2 states, but the same reasoning applies to 
other pairs of dipole-dipole coupled states. It is conve- 
nient to specify the atomic states by their orbital angu- 
lar momentum £ and azimuthal total angular momentum 
nij , i. e. as \£rrij), for example, |p + 3/2). The primary 
effect of the electric field is to lift the Zeeman degener- 
acy of the p 3 / 2 state so that W p ± 1 / 2 ^ W p ± 3 / 2 , although 



W pmj = W, 



p- mj . As a result, there are two sets of dipole- 



dipole coupled states with nearly degenerate R = oo 
asymptotes, which leads to the potential wells. The small 
spin-orbit splittings of the higher £ states should also lead 
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FIG. 1. (Color online) (a) The system under consideration 
consists of two Rydberg atoms. The relative position R of 
atom 2 with respect to atom 1 is expressed in terms of spher- 
ical coordinates. An external electric field E is applied in 
the z direction, (b) Internal level structure of each Rydberg 
atom. The Stark shift S = W p ±i/2 — W p ±- A /2 is positive for 
the above configuration. Only states connected by solid, blue 
dotted and red dashed lines are dipole coupled, (c) Potential 
curves for the M — 1 nsnp states as a function of the inter- 
nuclear spacing R. The Stark energy difference 5 between the 
\p ± 1/2) and \p ± 3/2) states is assumed to be positive. 



to wells, but their positions will be fixed. We define the 
energy splitting 8 by S = W p±1/2 - W p±3/2 . By the ap- 
propriate choice of AC or DC electric field 8 can be made 
positive or negative [13] • In Fig. [TJ; we show the 8 > 
case in which the \p ± 1/2) states lie above the \p ± 3/2) 
states at R = oo. 

The system of interest is the diatomic system of one 
nsx/2 atom and one np 3 / 2 atom, and we term its states 
the nsnp states. We describe the nsnp states as ordered 
direct products of atomic states. For example, atom 1 in 
the |s + 1/2) state and atom 2 in the \p+2>/2) state yields 
the state |s + l/2,p + 3/2). Constructing all the possible 
direct products of ns±/ 2 and np 3 / 2 states provides a com- 
plete set of basis functions for the nsnp states. With our 
choice of quantization axis the total azimuthal angular 
momentum M of the system remains a good quantum 
number in the presence of the electric field. There are 
four states of AI = 0, eight of M = ±1, and four of 
AI = ±2. We first consider the case in which the two 
atoms lie on the z axis, that is i2|[i3, and 8 > 0. The 
nsnp states have two asymptotic R = oo energies sepa- 
rated by 8, as shown by Fig. [TJ;. When R\\E M remains 
a good quantum number in the presence of the dipolc- 
dipole interaction, which is given by [l[ 

V dd = ±-[d (1) ■ dt 2) 3(d (1) • R)(d (2) ■ R)], (1) 



where d is the electric dipole-moment operator of atom 
i and R = R/R is the unit vector along the molecular 
axis. For any finite R the dipole-dipole interaction only 
couples states of the same M, and we focus on the four 
M = 1 states. We ignore the M = and ±2 states. 
There are four M = — 1 states degenerate with the AI = 1 
states, and the reasoning for the M = 1 states also ap- 
plies to them. As shown by Fig. [TJ;, as R — > oo there 
are two M = 1 states which converge to the \p ± 1/2) 
asymptote and two converging to the \p ± 3/2) asymp- 
tote. As R — > oo the former two states have dipole-dipole 
energy shifts proportional to ±l/i? 3 , while the latter pair 
are not dipole-dipole coupled to each other and have no 
first order dipole-dipole interaction. At small R, where 
Vdd 3> |<5| two of the four AI = 1 levels arc shifted up in 
energy and two down. Since the M = 1 levels can not 
cross, the lower level connected to the \p+l/2) asymptote 
must have a well, as shown in Fig. [JJ:, and the minimum 
occurs at a value of R such that Vdd ~ |<5|- For the eval- 
uation of Eq. (JTJ we compute the matrix elements of the 
dipole operator via the the Wigner-Eckart theorem [20[ 
according to 

v „ / 1- 1 

=v 

where C^2™, lq are Clebsch-Gordan coefficients and e q 
are orthonormal unit vectors arising from the decompo- 
sition of the dipole operator into its spherical compo- 
nents H^l- The reduced matrix element (np 3 / 2 \\d\\nsi/ 2 ) 
in Eq. ^ can be written in terms of a radial matrix el- 
ement between the np 3 / 2 and ns\/ 2 states 0, [2(| [2lj . 
We find V = \/l/3(np\r\ns), and for alkali atoms 
(np\r\ns) = n 2 for n ~ 40 [16|. Since the sum in Eq. ([2]) 
is a term of order unity, the magnitude of Vdd is given by 

il = \V\ 2 /R 3 . (3) 

Equating fi to \8\ yields the characteristic length Ro, 
given by 

Ro = (|2?|7l<5|) 1/3 = [n 4 /(3|£|)] 1/3 . (4) 

Fig. [JJ; also suggests that the depth of the potential well 
is approximately \8\/2. 

While it is clear that there is a potential well when 8 = 
0, two important questions remain. First, when 9^0 the 
dipole-dipole interaction couples states of different M, 
and we can no longer ignore states of M ^ 1 which cross 
the M = 1 states. What is the effect of these couplings? 
Second, what is the angular extent of the well in 61 Is 
it large enough to be useful? To address these questions 
we calculate the energy levels as functions of R and 9 by 
diagonalizing the Hamiltonian matrix resulting from the 
Stark shift 8 and the dipole-dipole interaction Vdd- We 
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FIG. 2. (Color online) (a)The 9 = M = ±1 potentials 
and the M = 0, ±2 potentials which cross them. The poten- 
tial well (black solid line) is two-fold degenerate. The blue 
dashed line corresponds to the M = ±2 states, and this en- 
ergy curve is two-fold degenerate. The red dashed-dotted line 
corresponds to a single M=0 state. The dotted green lines 
correspond to the van der Waals shifted M = ±1 states, and 
each of them is two-fold degenerate, (b) Expanded view of 
the point of five-fold degeneracy in Fig. [2^ for a small angle 
= 5°. Note that the figure looks qualitatively the same for 
even smaller angles. There is now a single state corresponding 
to the energy well (black solid line). 



first consider the effect of the other M states. The dipolc- 
dipole matrix elements have the following dependence on 
8. AM = 0: 1 - 3cos 2 0, AM = ±1: sin0cos0, AM = 
±2: sin 2 0. In general, all the levels are coupled. In 
Fig. [2k we show the calculated M = ±1 levels for = 
as well as the M = and ±2 levels which cross the 
well level. For clarity we have omitted the other M = 
and ±2 levels. In Fig. [2b we show the energy levels for 
8 = 5°. All the levels are coupled, and the degeneracies 
of the M = ±1 and M = ±2 states are broken. Equally 
clear, there are several potential wells, one of which has 
the same shape as the one shown in Figs. \Tjp and (2ji. 
Calculating the energy levels as a function of enables 
us to determine the angular extent of the well shown in 
Fig 12b, and in Fig. [3] we show a plot of the potential 
vs the x and z directions. The potential is azimuthally 
symmetric, and the potential wells exist along the +z and 
—z axes with widths of approximately 20°, wide enough 
to be useful. 

If we reverse the Stark shift, so that 6 < and the 
\p + 3/2) level lies above the \p + l/2) level, and consider 
8 = it/2, there is also a well. For 8 = tt/2, M = ±1 
states are coupled to each other but not to the M = 
and ±2 states. We again focus on the M = ±1 states, 
ignoring the M = and ±2 states. There are M = 
±1 states converging to both R — > oo asymptotes, and 
the well resembles the well of Figs. [TJ: and Ek, but it 
is nondegenerate in this case. Unlike the 9 = case, 
the strong AM = ±2 coupling lifts the degeneracy of 
the M = 1 and M = — 1 states. Due to the azimuthal 
symmetry the well is donut shaped as shown in Fig. [4] 
Calculating the energy levels as a function of 9 shows 
that the angular width of the potential about 9 = tt/2 is 
±20°, similar to the well shown in Fig. [3J 

Analytic diagonalization of the Hamiltonian matrix 



yields an expression for the potential wells shown in 
Figs.[2p andU Explicitly, 



1 



V(R) =-[35-412 
6 



(5) 



V280 2 + 35{n + 35 - 9n cos(20))] . 



This formula is valid for 5 > and S < 0. From the 
potential of Eq. ([S]) it is straightforward to derive expres- 
sions for the position of the minimum of the potential 
and the vibrational frequency of the radial motion in the 
potential. For = the minimum of the potential occurs 
at 



(T/e) 1 / 3 ^, 



(G) 



and the well depth is Ay = 0.57|<5|. The frequency of 
the vibrational motion is approximately given by 



Wvib 



(7) 



where (i is the reduced mass of the two atoms. 

For S < and = tt/2 the minimum of the potential 
occurs at 



= (4/3)^, 



(8) 



and the potential depth is AV = 0.75|<5|. The frequency 
of the vibrational motion is approximately given by 



w v ib = 2A\ 



Ij5\_ 
R 2 ^ 



(9) 



The Stark splitting of the level must be small com- 

pared to the np fine structure interval, which in Rb and 
Cs is ~ 0.01n~ 3 [22|. Accordingly, a reasonable limit for 
5 is \5\ < 1CT 4 ?i~ 3 . Combining this limit with Eq. (gj) 
yields Rq > (10 4 n/3) 1 / 3 n 2 , which is substantially larger 
than the Rydberg atom. 
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FIG. 3. (Color online) Potential of the 8 = well state in the 
x-z plane for 8 > 0. Two pronounced minima occur on the 
±z axes. The potential is azimuthally symmetric around the 
2 axis. 
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ergies just below the upper R = oo asymptote is purely 
attractive and can lead to collisional ionization; the other 
three have potential wells. Finally, it may be possible to 
detect the broadened time of flight field ionization signals 
due to atoms which remain close together 
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FIG. 4. (Color online) Potential of the 8 — -k/2 well state in 
the x-y plane for S < 0. The angular width of the well state 
around 8 — tt/2 is similar to the width of the corresponding 
well state for 8 > around 8 — 0, see Fig. [3] 



To develop a feeling for the dynamics in these shal- 
low potential wells, we use Rb n = 40 atoms with 
5/2ir = 10 MHz as an example. Combining these val- 
ues with Eqs. ©, @ and © leads to R = 8.25 x 10 4 , 
-Rmin = 8.69 x 10 4 , and w v ;b/27r = 24kHz, so motion on 
the potential curves is adiabatic. Since an equal linear 
superposition of the Rb 40s and 40p states has a K 
radiative lifetime of 120 /is, only three vibrational oscil- 
lations are likely [221 ]. 

The two atom states in these potential wells can be 
produced by laser excitation of pairs of atoms to the 
diatomic nsns state with subsequent excitation to the 
dipole-dipole bound nsnp state by microwave excitation. 
While the most interesting prospect is forming the dipole- 
dipole wells to match the spacing of atoms in a lat- 
tice [I?} , it is useful to consider other ways in which these 
wells can be manifested. Definitive evidence can be ob- 
tained by taking advantage of the fact that in zero field 
the potentials are either attractive or repulsive. There 
are no long range wells in which atoms can be trapped, 
and atoms are not stable at energies slightly removed 
from the R = oo energies. Stable atoms in these poten- 
tial wells can be detected in several ways. First, driving 
a second microwave transition from the nsnp well to, for 
example, the nsn'd state, which has no dipole-dipole in- 
teraction, should show a satellite feature removed from 
the atomic transition by the well depth. This feature 
should become more pronounced if a time delay is intro- 
duced before the second microwave transition, the lim- 
iting factor being radiative decay. Second, the motion 
of atoms on attractive dipole-dipole potentials in zero 
field leads to collisional ionization on time scales much 
faster than the radiative lifetime (23l - [26j . This ionization 
should be significantly suppressed by the existence of the 
long range potential dipole-dipole wells, for example, in 
the = tt/2 case with S < 0, which leads to the donut 
potential of Fig. [H only one of the four potentials at en- 



In conclusion, the Stark tuned dipole-dipole potentials 
described here provide a novel, adjustable way to locate 
interacting atoms at a chosen distance. 
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